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;_ i ■ Abstract: In this paper, we propose a reproducing kernel Hilbert space 

Q-f method (RKHSM) for solving the sine-Gordon (SG) equation with initial and 

^ ■ boundary conditions based on the reproducing kernel theory. Its exact solution 

p^ \ is represented in the form of series in the reproducing kernel Hilbert space. 

Some numerical examples have been studied to demonstrate the accuracy of the 
present method. The results obtained from the method are compared with the 
exact solutions and the earlier works. Results of numerical examples show that 
the presented method is simple and effective. 
(-h Keywords: Reproducing kernel method, series solutions, sine-Gordon equa- 

■^ ' tion, reproducing kernel space. 
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1. Introduction 

Since Russell's [1-2] first observation for some special waves with charac- 
teristic properties on a canal in 1834 many scientists among whom Korteweg 
£> | and de Vries [3] , who derived the equation concerning the propagation of waves 

in one direction on the free surface of a shallow canal also known as the KdV 
equation, have investigated them more extensively. The name soliton was first 
used in Zabusky and Kruskal [4] in order to emphasize that a soliton is a lo- 
calized entity which keeps its identity after interaction. Equations which also 
lead to solitary waves are the sine-Gordon, the cubic Schrodinger equation, etc. 
Because of the great importance of the study of physical phenomena theoret- 
ical solutions of nonlinear and especially of soliton type equations have been 
developed during the last years [5]. 
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The nonlinear one-dimensional Sine-Gordon (SG) equation is a very im- 
portant nonlinear hyperbolic partial differential equation PDE. It appears in 
d ■ differential geometry gained its significance because of the collisional behaviors 

of solitons that arise from these equations. It was basically considered in the 
nineteenth century in the course of surface of constant negative curvature. This 
equation attracted a lot of attention in the 1970s due to the presence of soliton 
solutions [6-7]. The sine-Gordon equation comes out in a number of physical 
applications [8-10] including applications in the chain of coupled pendulums 
and modelling the propagation of transverse electromagnetic (TEM) wave on a 
superconductor transmission system. 

Consider the one-dimensional nonlinear sine-Gordon equation 



dt 2 ""' " dx 



"a^fo*) = ol^fo*) ~ sin(u(x,t)), a<x<b,t>0, (1.1) 



with initial conditions 

u(x,0) = f(x), a<x<b, (1.2) 

du 

— (x,Q) = g{x), a<x<b, 

and the boundary conditions 

u(a,t) = hi(t), u(b,t) = h 2 {t), t > 0. (1.3) 

Since the exact solution of the SG equation can only be obtained in special 
situations, many numerical schemes are constructed to solve the SG equation. 
The study of numerical solutions of the sine-Gordon equation have been inves- 
tigated considerably in the last few years. For solving (1.1), for instance, high- 
order solution of one-dimensional sine-Gordon equation using compact finite 
difference and DIRKN methods [8] . The authors of [11] proposed a numerical 
method form solving (1.1) using collocation and radial basis functions. Also, the 
boundary integral equation approach is used in [12] . Bratsos has proposed a 
numerical scheme for solving one dimensional sine-Gordon equation and a third- 
order numerical scheme for the two-dimensional sine-Gordon equation in [13-14] 
respectively. Also, in [15] , a numerical method using radial basis function for the 
solution of two-dimensional sine-Gordon equation is used. In addition, several 
authors recommended spectral methods and Fourier pseudospectral method for 
solving nonlinear wave equation using a discrete Fourier series and Chebyshcv 
orthogonal polynomials [16-18]. Ma and Wu [19] presented a meshless scheme 
by using a multiquadric (MQ) quasi-interpolation named L £ without solving 
a large-scale linear system of equations, but a polynomial p(x) was needed to 
improved the accuracy of the scheme. More recently, another numerical work 
has been investigated by Khaliq et al. [20] . 

One difficult point of numerically solving the sine-Gordon equation is how 
to solve the nonlinear system resulting by discretization. The Newton iteration 
method is often used as a basis for designing numerical schemes. In this case, the 
Jacobian matrices have to be established, inverted and possibly updated during 
the iteration. When high approximation accuracy is desired, it requires one to 
use a sufficiently small time-step and sufficiently fine grids. Thus it demands a 
large amount of computational effort [21]. 

In this paper, we solve Eqs. (1.1) and (1.3) by using Reproducing Kernel 
Method (or RKM). The nonlinear problem is solved easily and elegantly by using 
RKM. The technique has many advantages over the classical techniques. It 
also avoids discretization and provides an efficient numerical solution with high 
accuracy, minimal calculation, avoidance of physically unrealistic assumptions. 
In the next section, we will describe the procedure. 

The theory of reproducing kernels [22] was used for the first time at the be- 
ginning of the 20 th century by Zaremba in his work on boundary value problems 



for harmonic and biharmonic functions. Reproducing kernel theory has impor- 
tant applications in numerical analysis, differential equations, probability and 
statistics and so on [23-35]. Recently, using the RKM, some authors discussed 
fractional differential equation, nonlinear oscillator with discontinuity, singular 
nonlinear two-point periodic boundary value problems, integral equations and 
nonlinear partial differential equations and so on [23-35]. 

The paper is organized as follows. Section 2 is devoted to several reproducing 
kernel spaces. Solution represantation in W (O) and a linear operator have 
been presented in Section 3. Section 4 provides the main results, the exact 
and approximate solution of Eqs.(l.l) and (1.3) and an iterative method are 
developed for the kind of problems in the reproducing kernel space. We have 
proved that the approximate solution converges to the exact solution uniformly. 
Some numerical experiments are illustrated in Section 5. We provide some 
conclusions in the last section. 

2. Preliminaries 

2.1. Reproducing Kernel Spaces 

In this section, we define some useful reproducing kernel spaces. 

Definition 2.1. (Reproducing kernel). Let E be a nonempty abstract set. 
A function K : E x E — > C is a reproducing kernel of the Hilbert space H if 
and only if 

a) Vt€E, K(.,t) eH, 

b)VteE,^EH, (<p(.),K(.,t))=<p(t). 

The last condition is called "the reproducing property" the value of the 
function if at the point t is reproduced by the inner product of tp with K (.,£). 

Definition 2.2. 

{u(x) | u(x), u'(x), u"{x) are absolutely continuous real value functions in [0, 1], 
u^{x) G L 2 [0, 1], x G [0, 1], w(0) = 0, w(l) = 0, 
The inner product and the norm in W$ [0, 1] are defined respectively by 

,4 

(u(x),g(x)) w3 = u(0)g(0)+u'(0)g'(0)+u'(l)g'(l)+ / u^(x)g^(x)dx, u(x), g(x) G W 2 3 [0,1], 

Jo 

and 



h\\wi = yJM wS , «ew![o,i] 



The space W% [0, 1] is a reproducing kernel space, i.e., for each fixed y 6 [0, 1] 
and any u(x) € Wf [0, 1], there exists a function R y (x) such that 

u(y) = (u(x), R y (x)) W 3 . 

Definition 2.3. 

(v(t) | u(t), u'(t),u"(t) are absolutely continuous in [0,T], 
u^(t) Gi 2 [0,l], te [0,T], w(0) =0, u'(0) = 0. 
The inner product and the norm in W^[0,T] are defined respectively by 

(u(t),g(t)) 8 =2«W(0)flW( )+ / u< 3 \t)g^(t)dt, u(t),g(t) € W 2 3 [0,T], 

and 



||«|| w , = JW^~ r ueWi[0,T]. 



The space W 7 ^ [0,T] is a reproducing kernel space and its reproducing kernel 
function r s (t) is given by 



1 „2+2 i 1 „2^3 L, f 4j 1 ^5 

rs (t) = 



• O ij I i o 04. I - 1 on ) — ) 



. 6 t- i 1 9 ^ v 24 ~ 120 ' 



Definition 2.4. 

H/1 r n 11 _ / ^0*0 I w(#) i s absolutely continuous in [0, 1], 

2[ ' J ~\ u'{x) ei 2 [0,l], xe [0,1]. 

The inner product and the norm in W\ [0, 1] are defined respectively by 

(u(x),g(x)) wl =u(0)g(0)+ I u'(x)g'(x)dx, u(x),g(x) € W%[0, 1] 

Jo 

and 



yjM wi > uew^o.i]. 



The space Wi^O, 1] is a reproducing kernel space and its reproducing kernel 
function Q y (x) is given by 

{1 + x, x < y, 
l + y, x>y. 

Definition 2.5. 



{u(i) | u(t) is absolutely continuous in [0, T], 
u'(t)eL 2 [0 7 T], te[0,T\. 
The inner product and the norm in W^O, T] are defined respectively by 

(u(t),g(t)) = u(0)g(0)+ [ u'(t)g'(t)dt, u(t),g(t) € W}[0,T], 

W 2 JO 

and 



NlwJ = \J( U ' U )w^ U ^ W 2\®> T ] 



The space W\ [0, T] is a reproducing kernel space and its reproducing kernel 
function q s (t) is given by 

( 1 + t, t<s, 
I 1 + s, t > s. 



Theorem 2.1. The space W% [0, 1] is a complete reproducing kernel space 
whose reproducing kernel function R y (x) is given as, 

EiLiCiG/)^- 1 , x < y, 

H,U)={ ( (2.1) 



where 



ci(y) = o, 

1 5 5 4 127 2 31 

^ vyy 122 y 244 y 244 y 61 y 

1 5 127 4 1 , 1137 2 127 

c 3 (v) = y° H y y 6 -\ y y, 

Ky ' 2928 y 5856 y 12 y 1952 ,y 244 y ' 

c 4 (y) = 0, 

c 5(yj - 2 938 y 5856 y 5856 17 1464^' 

ceU/j - 732Q y + 2928 y 2928 ?y 122 y+ 120' 

dl(y) ^i4o y5 ' 

J . . 1 5 31 4 127 2 31 

vyy 122 ,y 1464 y 244 y 61 y 

1 5 127 4 1137 2 127 

yy ' 2928 y 5856 y 1952 y 244 y ' 
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1 5 5 4 127 2 5 

dJv) = y b y 4 H y z H y, 

byyj 2928 tf 5856 w 5856 y 244 y ' 

1_ 5 _J_ 4 1 2 1 

eu/j - 732Q y + 2g28 2/ 2g28 y 122 2/- 

Proof: Let u 6 W 2 3 [0, 1] and < y < 1. Define R y by (2.1). Note that 

f ELi^+iO/) 2 ^ 1 ; x < y> 
i?;(x) = <^ 

[ ELi*(* + 1 ) c '+2(y) x ^ 1 . x < v-. 

I ELi*(* + 1 )^+2 (y)^" 1 , x > y, 
{ Eti^ + ^ + ^+^y)^- 1 , x<y, 

{ Eti i(* + !)(* + 2)ci 4+3 (y)x 4 - 1 , x > j/, 
f ELi J ( J + 1 )( J + 2 )(* + 3)c t +4(2/)^- 1 , x<y, 

4 4) (*) = 

I ELi *(* + !)(* + 2)(z + 3)d i+4 (j/)*- 1 , x > y, 

and 

C 120c 6 (y), x < y, 

4 5) (*) = { 

{ 120d 6 (y), x>y. 
By Definition 2.2 and integrating by parts two times, we obtain 

(u(x),R y (x)) wi = u(0)Ry(0) + u'(0)R' y (0) + u'(l)R' y (l) + u"(l)R ( y 3 \l) 

2 +u(l)R y 5 \l)-u(0)R y 5 \Q)-tiu(x)R y 6 \x)dx 



= u'(0)(R' y (0) + 4 4) (0)) + u'(l)(R' y (l) - R y 4 \l)) 
+u"(l)R y 3) (l)-u"(0)R y 3 \0) 
+ Jo R y 5) (x)u(x)dx+ j 1 R y 5) (x)u(x)dx 

= «'(0)(ca(i/) + 24c 5 (y)) - u"(0)(6c 4 (y)) 

+ u'(l){d 2 (y) + 2d 3 (y) + MM - 20d 5 (j/) - 115d 6 (y)) 

+«"(l)(6d4(|/) + 2Ad 5 {y) + 60d 6 (y)) 
+ f£ 120c 6 (y)u(x)dx+ J 1 120d 6 (y) u (x)dx 

= 120u(y)(^) = u(y). 



Definition 2.6. 



W(fi): 



u(x,t) I a^dt* ' ^ S completely continuos in O = [0, 1] x [0,T], 



J^ e L 2 (ft) , u(x, 0) = 0, ^1 = o, „( 0) t) = 0, «(1, t) = 0. 
The inner product and the norm in W (ft) are defined respectively by 



2 - T 



(u(x,t),g(x,t))^ 



E 



/=o 



a 3 a* 



a 3 a* 



^ 3 ^ u(0 ' <) ^a^ 5(0 ' t) 



r// 



*-- / a-? a? 



j=0 



dP 



dp- 



W: 



I /-T r 9 3 £3 £3 ^3 



o Jo 



dtdx, 



and 



Nl w = vW>w" we VK(ri) 



Theorem 2.2. W(ft) is a reproducing kernel space and its reproducing 
kernel function is 

K (y,s) ( x ,t) = R y (x)r s (t), 

such that for any u(x, t) £ W (ft) , 



and 



u(y, s) = (u(x, t),K(y >s ) (x, t)) w , 

K {y.s) (X,t) =K(x,t) (V, s )- 



Definition 2.7. 

^~ C u(x,t) | u(x,t) is completely continuos in fl= [0, 1] x [0,T], 

The inner product and the norm in W (ft) are defined respectively by 



(u(x,t),g(x,t)) 



w 



T v 



m uM m gM 



+ {u(x,0),g(x,0) w 

,1 rT r g q 



(II 



n, at u{Xit) mm g{Xit) 



dtdx. 



and 



Ml 



w 



(u,u)^, ueW(Cl). 



W (CI) is a reproducing kernel space and its reproducing kernel function Gi ys \ (x, t) 
is 

G(v,a) (#,*) = Qy(x)q s (t). 



3. Solution represantation in W (CI) . 

In this section, the solution of equation (1.1) is given in the reproducing 
kernel space W (CI). On defining the linear operator L : W (CI) — > W (CI) as 

d 2 v . . d 2 v . . 

After homogenizing the initial and boundary conditions model problem (1.1) 
changes the following problem: 



Lv = M(x, t, v(x, t)), (x, t) e [0, 1] x [0, T], 

v(x,0) = §£(x,0) = v(0,t) = v(l,t) = 0. 
We replace v(x,t) with u(x,t) in (3.1), for simplicity. 

Lemma 3.1. L is a bounded linear operator. 

Proof: 

,Tr a 



(3.1) 



Lu\ 



w 



dt Lu(0,t) 

1 C f aa 

y [a^at 



|^(o, *) 



dt + (Lu(x, 0), Lu(x, 0)) w i 



Lu(x, t) 



dxdt 



dt+[Lu(Q,0)Y 



1 2 



— Lu(x,0) 



dx 



T 

o Jo 



Txm LuM 



dxdt, 



since 



u(x,t) = (u(£,ri),K {Xit) (£,ri)) w , 
Lu(x,t) = (u(£,ri),LK( Xtt) (t,ri)) w , 

from the continuity of Ki xt \ (£, 77), we have 

\Lu(x,t)\ < \\u\\ w \\LK( x>t) (tv)\\ w < «o \\u\\w- 



similarly for i = 0, 1 



Qi I Qi 

^—Lu{x,t) = { u(£,»7), ^iK(i )t ) (£,17) 



dx 

d_d_ 

dt dx 1 



dx 1 
d d l 



\Y 



Lu{x, t) = ( «(£, rf), ——-LK {Xit) (£, rf) 



and then 



— — rLu(x, t) 

dx 1 v ' ' 



dt dx 1 



<ei\\u\\ w . 



\Y 



d d i 
dt dx' 



Lu{x, t) 



</ill«ll W . 



Therefore 



\Lu(x,t)\\h <E( e i+tf)N&<« 2 NI 



w ■ 



a 



i=0 



Now, choose a countable dense subset {(x\,t\) , (X2, £2) , •••} in O = [0, 1] x 
[0, T] and define 

<Pi(x, t) = G( Xiiti ) (x, t) , $>i(x, t) = L*<p t (x, t), 

where L* is the adjoint operator of L. The orthonormal system < $>i(x,t) > 

of W (f2) can be derived from the process of Gram-Schmidt orthogonalization 
oi {*i(x,t)}Zi & 

i 

$ 4 (x,i)=^/3 lfc * fc (x,i). 
fc=i 



Theorem 3.1. Suppose that {(xi,ti)} c *l 1 is dense in il; then {^"i(x, i)}^i 
is complete system in W (ft) and 

*i(a;,t) - L {VtB) K M (aJ,*)| (l/ia)=(a . liti) • 

Proof: We have 

$i(x,t) = {L*<p i )(x,t) = ((L*v i )(y,s),K {Xit) (y,s)) w 
= (<Pi(y,s), L( y>3) K( Xtt) (y,s))ft 

= L (y,s)K {x ,t){y,s)\ {ys)={XtU) 

Clearly * 4 (x, i) G VF (O) . For each fixed u(x, i) <E VK (ft) , if 
(«(a;,t),*i(a;,t)> W r = ) i = 1,2,... 



then 

(u{x,t),(L*<p i )(x,t)) w = (Lu(x,i),<Pi(x,i))ft 

= (Lu) (xi, U) = 0, z = 1,2,.... 

Note that {(xi, U)}°1 1 is dense in ft, hence, (Lu) (x, t) = 0. It follows that u = 
from the existence of L . So the proof is complete. □ 

Theorem 3.2. If {(xi,£i)}Si is dense in £7, then the solution of (3.1) is 

oo i 

u{x,t)=J2J2 /3 ik M (xk, t k , u(x k , t k ))h(x,t). (3.2) 

i=l fc=l 

Proof. Since {^,(x, t)}°l 1 is complete system in W (ft) , we have 

oo 
U(x,t) = y2(u(x,t),$i(x,tj) $i(x,t) 

, — , \ f w 

oo i 
= J2J2^ik( U ( X ^)^k(x,t)) w ^i(x,t) 

i=\ fe=l 

oo % 

= X) X)^ M*' *)' L *Vk(x, *)) W ^i( X > f ) 
i=l fe=l 

oo i 

= ^2^20ik( L u(x,t),ip k (x,t))^^i(x,t) 
i=\ fc=l 

oo z 

i=i fc=i 

oo i 

= EE^ 1 *' 1 ' 1 ' 1 )*^*) 
i=i fc=i 

oo i 

= EE^f^' ife ' u ( xk > ^fe))**^)*)- 
i=i fc=i 

Now the approximate solution u n (x,t) can be obtained from the n— term 
intercept of the exact solution u(x, t) and 

u n (x,t) ^^^fcMfn, i fe , u(a;fe, *&))** (a;,*). 



i=i fc=i 
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Obviously 

\\u n (x,t) — u{x, t)\\ — *• 0, (n -4 oo) . D 

4. The method implementation 

If we write 

i 

A=^2(3 lk M(x k , t k , u(x k , t k )), 

k=l 

then (3.2) can be written as 

oo 

u(x,t) =^ j A i $ i (x,t). 

Now let (xi,ti) = 0; then from the initial and boundary conditions of (3.1), 
u(x,\,t\) is known. We put uo(xi,ii) = u(x\.,t\) and define the n— term 
approximation to u(x, t) by 

n 

«„(x,t) = 5^Bi$i(a;,t). (3.3) 

where 

Bi = ^f3 ik M(x k , t k , u k -i(x k , tk))- (3.4) 

/c=i 

In the sequel, we verify that the approximate solution u n (x,t) converges to 
the exact solution, uniformly. 

Theorem 4.1. Suppose that ||u n || is a bounded in (3.3) and (3.1) has a 
uniqe solution. If {(2^,^)}°^ is dense in W (Q) , then the n— term approxi- 
mate solution u n {x,t) can be derived from the above method converges to the 
analytical solution u(x,t) of (3.1) and 

00 
u(x,t) =^B,$ 4 (x,i), 

i=l 

where Bi is given by (3.4). 

Proof: First, we prove the convergence of u n {x, t). From (3.3), we infer that 
u n+ i(x,t) = u n (x,t) + B Il+1 * n+ i(x,i), 

The orthonormality of l^if yields that 

n+1 

\\u n+l \\ 2 = \\u n \\ 2 +B 2 n+1 =Y J BJ. (3.5) 



i=\ 



11 



In terms of (3.5), it holds that ||u„ + i|| > \\u n \\ . Due to the condition that ||w„ 
is bounded, ||u„|| is convergent and there exists a constant c such that 

oo 

This implies that 

{Bi}Zi e i 2 - 

If 77i > n, then 

||«m-Wn|| = \\u m -Um-l+U m -l—Um-2 + ---+Un+l-U n \\ 

< \\u m -u m -i\\ + \\u m -i - u m - 2 \\ + ... + ||Wn+i - Wn|| • 
On account of 

||Wm-Um-l|| = B 2 m , 

consequently, 



i n \\ = 2_. B? — > 0, as n — >• oo. 



| ^m ^n , , 

l=n+l 

The completeness of W (fl) shows that u n — > u as n -4 oo . Now, let we prove 
that u is the solution of (3.1). Taking limits in (3.3) we get 

oo 

u(x,t) = y^ B$i(x,t). 

i=l 

Note that 

OO 

(Lu)(x,i) =^ £;!,$; (Z,i), 
i=l 

and 

oo 

(Lu)(xi,ti) = ^2BiL^i(xi,ti) 



IV 



vv 



2—1 

Therefore 

z oo / i \ 

Y / f3 a (Lu)(x l ,t l ) = J2Bi(*i(x,t), ^0 a *i(x,t) 
i=i »=i \ ;=i / 

oo 



oo 

J2Bi($i(x,t), L*^(x,t) 

i=l 

oo 

J2Bi($i(x,t), *i(x,t) 



i—1 
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In wiev of (3.4), we have 

Lu(x l ,t l ) = M{x v ti, u(xi, ti)) 

Since {{xi,ti)}°l 1 is dense in $7, for each (y, s) £ 0, there exists a subsequence 
{ (xnj , t nj ) }°l x such that 

(x nj ,t ni ) -> (y,s), (j -> oo). 
Wc know that 

Lu\x njl t n .j =Myx n - 1 t njl u(x nj , t n .)). 

Let j — > co; by the continuity of /, we have 

(Lu)(y, s) = M(y, s, u(y, s)). 

which indicates that u(x, t) satisfies (3.1). 

Remark 4.1. In a same manner, it can be proved that 
du n (x,t) du(x,t) 



— >• 0, as n — > oo, 



where 



and 



dx dx 

du{x,t) ^ D 9*,(x,i) 



dx 

du n (x,t) 
dx 



i=i 






dx 



Bi is given by (3.4). 



5. Numerical Results 

Example 5.1. In this example, we consider SG equation (1.1) without 
nonlinear term sin(u) in the region [0,1] x [0,T]. The initial conditions are 
given by 

du 
u(x, 0) = sin(7ra;), — -(x, 0) = 0, 



at 



with the boundary conditions 



u(0,t) =u(l,t) = 0. 
The exact solution is given in [36] as 

U (x, t) — — (sin(7r(x + t) + sin7r(x — t)). 
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After homogenizing the initial and boundary conditions we obtain (5.1) as 

(5.1) 



°£(x,t)-%$(x,t) = -n 2 sm{Trx), 



u(x, 0) = §£(a;, 0) = u(0, t) = u(l,t) = 0. 

Example 5.2. In this example, we take notice of SG equation (1.1) in the 
region [0, 1] x [0, T]. The initial conditions are given by 

du 



w(a;, 0)=0, — (x, 0) = 4secft(x), 



The exact solution is 



u{x,t) — arctan(sec h(x)t). 
After homogenizing the initial conditions we obtain (5.2) as 



cYu 
~dt T 



41 , 8£ sinh (x) 



(x,t) &{x,t) =- S m(u(x,t)+Usech(x)) - ^ + ^^ 
u(x, 0) = ||(x, 0) = u(0, t) = u(l,t) = 0. 



(5.2) 



X 


t 


Exact Solution 


Approximate Solution 


Absolute Error 


Relative Error 


Time 


0.1 


0.1 


0.2938926262 


0.2938930965 


4.703 x 10-' 


1.6002443 x 10- b 


3.860 


0.2 


0.2 


0.4755282582 


0.4755313577 


3.0995 x 10- b 


6.518014327 x 10~ b 


3.016 


0.3 


0.3 


0.4755282582 


0.4755183355 


9.9227 x 10- b 


2.086668842 x 10~ b 


2.984 


0.4 


0.4 


0.2938926261 


0.2939007109 


8.0848 x 10~ 6 


2.750936663 x 10~ b 


3.000 


0.5 


0.5 


0.0 


0.0000282140 


2.82140 x 10~ 5 


oo 


3.094 


0.6 


0.6 


-0.2938926264 


-0.2939063137 


1.36873 x 10- b 


4.657245120 x 10~ 5 


3.031 


0.7 


0.7 


-0.4755282583 


-0.4755305759 


2.3176 x 10~ 6 


4.8737377 x 10~ b 


3.031 


0.8 


0.8 


-0.4755282581 


-0.4755277748 


4.833 x 10-'' 


1.016343386 x 10~ b 


2.953 


0.9 


0.9 


-0.2938926260 


-0.2938966580 


4.0320 x 10~ 6 


1.371929624 x lO" 5 


3.204 


1.0 


1.0 


0.0 


-3.690702068 x 10" Y 


3.690702068 x 10~' 


oo 


3.578 



Table 1. Numerical solutions for Example 5.1. 



X 


ES 


AS 


AE 


RE 


Time CPU(s) 


-0.80 


0.5877852522 


0.5877854278 


1.756 x 10~ 7 


2.987485639 x 10" Y 


1.404 


-0.40 


0.9510565165 


0.9510565935 


7.70 x 10-* 


8.096259125 x 10~ 8 


1.373 


0.40 


-0.9510565165 


-0.9510565935 


7.70 x 10-* 


8.096259125 x 10~ 8 


1.373 


0.80 


-0.5877852522 


-0.5877854278 


1.756 x 10~ 7 


2.987485639 x 10" Y 


1.404 
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Table 2. Numerical solutions for example 5.1 for t = 1. 



X 


AE [37] 


AE [RKHSM] 


RE [37] 


RE [RKHSM] 


-0.80 


1.94E-05 


1.756£-07 


3.29£-05 


2.987485639^-7 


-0.40 


2.84£ - 07 


7.700£ - 08 


2.98£-07 


8.096259125^-8 


0.40 


2.84£ - 07 


7.700£ - 08 


2.98£-07 


8.096259125^-8 


0.80 


1.94E-05 


1.756£-07 


3.29£-05 


2.987485639^-7 



Table 3. Comparison Absolute error and Relative error for Example 5.1. 



X 


t 


Exact Solution 


Approximate Solution 


Absolute Error 


Relative Error 


Time 


0.1 


0.1 


0.09917563307 


0.09917562 


1.307 x 10~ M 


1.317864035 x 10" Y 


1.688 


0.2 


0.2 


0.1936096360 


0.1936060475 


3.5885 x 10~ b 


1.853471797 x 10~ b 


1.719 


0.3 


0.3 


0.2794771925 


0.2794797547 


2.5622 x 10~ b 


9.167832184 x 10~ b 


0.593 


0.4 


0.4 


0.3543825410 


0.354382187 


3.540 x 10~ 7 


9.989205422 x 10~ 7 


0.608 


0.5 


0.5 


0.4173597173 


0.417357842 


1.8753 x 10~ b 


4.493246287 x 10~ b 


0.577 


0.6 


0.6 


0.4685399032 


0.468538681 


1.2222 x 10~ b 


2.60852916 x 10~ b 


0.577 


0.7 


0.7 


0.5087309794 


0.508730850 


1.294 x 10~ 7 


2.543584040 x 10~ 7 


0.593 


0.8 


0.8 


0.5390654126 


0.539066210 


7.974 x 10"' 


1.47922679 x 10~ b 


0.608 


0.9 


0.9 


0.5607645949 


0.560770310 


5.7151 x 10~ b 


1.019162061 x lO" 15 


0.546 


1.0 


1.0 


0.5750061826 


0.575005129 


1.0536 x 10~ b 


1.832328124 x 10~ b 


0.561 



Table 4. Numerical solutions for Example 5.2. 



X 


ES 


AS 


AE 


RE 


Time CPU(s) 


-0.80 


2.568109722 


2.568109726 


4 x 10- 9 


1.557565849 x 10~ 9 


0.733 


-0.40 


2.985843344 


2.985843344 


0.0 


0.0 


0.686 


0.00 


3.141592654 


3.141592654 


0.0 


0.0 


0.702 


0.40 


2.985843344 


2.985843344 


0.0 


0.0 


0.686 


0.80 


2.568109722 


2.568109726 


4 x 10~ 9 


1.557565849 x 10~ 9 


0.733 



Table 5. Numerical solutions for example 5.2 for t = 1. 



X 


AE [37] 


AE [RKHSM] 


RE [37] 


RE [RKHSM] 


-0.80 


1.53.E-08 


4£-09 


5.96£ - 09 


1.557565849S- 


-09 


-0.40 


3.54S-10 


0.0 


1.18.E-10 


0.0 


0.00 


1.62S-10 


0.0 


5.15^-11 


0.0 


0.40 


3.54S-10 


0.0 


1.18.E-10 


0.0 


0.80 


1.53S-08 


4£-09 


5.96£ - 09 


1.557565849S- 


-09 



Table 6. Comparison Absolute error and Relative error for example 5.2. 
6. Conclusion 



15 



In this study, linear and nonlinear SG Equations were solved by RKHSM. 
We described the method and used it in some test examples in order to show 
its applicability and validity in comparison with exact and other numerical so- 
lutions. The obtained results show that this approach can solve the problem 
effectively and need few computations. The results are satisfactory. The results 
that we obtained were compared with the results that were obtained by [37]. 
Numerical experiments on test examples show that our proposed schemes are of 
high accuracy, and support the theoretical results. According to these results, 
it is possible to apply RKHSM to linear and nonlinear differential equations 
with initial and boundary conditions. It has been shown that the obtained 
results arc uniform convergent and the operator that was used is a bounded lin- 
ear operator. Thus the studies that were compared with this method and this 
study show that RKHSM can be apply to high dimensional partial differential 
equations, integral equations and fractional differential equations without any 
transformation or discretization and good results can be obtained. 
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